A novel variant in GLIS3 is associated with osteoarthritis

Objectives Osteoarthritis (OA) is a complex disease, but its genetic aetiology remains poorly characterised. To identify novel susceptibility loci for OA, we carried out a genome-wide association study (GWAS) in individuals from the largest UK-based OA collections to date. Methods We carried out a discovery GWAS in 5414 OA individuals with knee and/or hip total joint replacement (TJR) and 9939 population-based controls. We followed-up prioritised variants in OA subjects from the interim release of the UK Biobank resource (up to 12 658 cases and 50 898 controls) and our lead finding in operated OA subjects from the full release of UK Biobank (17 894 cases and 89 470 controls). We investigated its functional implications in methylation, gene expression and proteomics data in primary chondrocytes from 12 pairs of intact and degraded cartilage samples from patients undergoing TJR. Results We detect a genome-wide significant association at rs10116772 with TJR (P=3.7×10−8; for allele A: OR (95% CI) 0.97 (0.96 to 0.98)), an intronic variant in GLIS3, which is expressed in cartilage. Variants in strong correlation with rs10116772 have been associated with elevated plasma glucose levels and diabetes. Conclusions We identify a novel susceptibility locus for OA that has been previously implicated in diabetes and glycaemic traits

The 2018 Otto Aufranc Award: How does genome-wide variation affect osteolysis risk after THA?

BACKGROUND: Periprosthetic osteolysis resulting in aseptic loosening is a leading cause of THA revision. Individuals vary in their susceptibility to osteolysis and heritable factors may contribute to this variation. However, the overall contribution that such variation makes to osteolysis risk is unknown. QUESTIONS/PURPOSES: We conducted two genome-wide association studies to (1) identify genetic risk loci associated with susceptibility to osteolysis; and (2) identify genetic risk loci associated with time to prosthesis revision for osteolysis. METHODS: The Norway cohort comprised 2624 patients after THA recruited from the Norwegian Arthroplasty Registry, of whom 779 had undergone revision surgery for osteolysis. The UK cohort included 890 patients previously recruited from hospitals in the north of England, 317 who either had radiographic evidence of and/or had undergone revision surgery for osteolysis. All participants had received a fully cemented or hybrid THA using a small-diameter metal or ceramic-on-conventional polyethylene bearing. Osteolysis susceptibility case-control analyses and quantitative trait analyses for time to prosthesis revision (a proxy measure of the speed of osteolysis onset) in those patients with osteolysis were undertaken in each cohort separately after genome-wide genotyping. Finally, a meta-analysis of the two independent cohort association analysis results was undertaken. RESULTS: Genome-wide association analysis identified four independent suggestive genetic signals for osteolysis case-control status in the Norwegian cohort and 11 in the UK cohort (p ≤ 5 x 10). After meta-analysis, five independent genetic signals showed a suggestive association with osteolysis case-control status at p ≤ 5 x 10 with the strongest comprising 18 correlated variants on chromosome 7 (lead signal rs850092, p = 1.13 x 10). Genome-wide quantitative trait analysis in cases only showed a total of five and nine independent genetic signals for time to revision at p ≤ 5 x 10, respectively. After meta-analysis, 11 independent genetic signals showed suggestive evidence of an association with time to revision at p ≤ 5 x 10 with the largest association block comprising 174 correlated variants in chromosome 15 (lead signal rs10507055, p = 1.40 x 10). CONCLUSIONS: We explored the heritable biology of osteolysis at the whole genome level and identify several genetic loci that associate with susceptibility to osteolysis or with premature revision surgery. However, further studies are required to determine a causal association between the identified signals and osteolysis and their functional role in the disease. CLINICAL RELEVANCE: The identification of novel genetic risk loci for osteolysis enables new investigative avenues for clinical biomarker discovery and therapeutic intervention in this disease

Allelic expression analysis of the osteoarthritis susceptibility locus that maps to MICAL3

<p>Abstract</p> <p>Background</p> <p>A genome-wide association scan with subsequent replication study that involved over 67,000 individuals of European ancestry has produced evidence of association of single nucleotide polymorphism rs2277831 to primary osteoarthritis (OA) with a P-value of 2.9 × 10^-5. rs2277831, an A/G transition, is located in an intron of <it>MICAL3</it>. This gene is located on chromosome 22q11.21 and the association signal encompasses two additional genes, <it>BCL2L13 </it>and <it>BID</it>. It is becoming increasingly apparent that many common complex traits are mediated by <it>cis</it>-acting regulatory polymorphisms that influence, in a tissue-specific manner, gene expression or transcript stability.</p> <p>Methods</p> <p>We used total and allelic expression analysis to assess whether the OA association to rs2277831 is mediated by an influence on MICAL3, BCL2L13 or BID expression. Using RNA extracted from joint tissues of 60 patients who had undergone elective joint replacement surgery, we assessed whether rs2277831 correlated with allelic expression of either of the three genes by: 1) measuring the expression of each gene by quantitative PCR and then stratifying the data by genotype at rs2277831 and 2) accurately discriminating and quantifying the mRNA synthesised from the alleles of OA patients using allelic-quantitative PCR.</p> <p>Results</p> <p>We found no evidence for a correlation between gene expression and genotype at rs2277831, with P-values of 0.09 for <it>BCL2L13</it>, 0.07 for <it>BID </it>and 0.33 for <it>MICAL3</it>. In the allelic expression analysis we observed several examples of significant (p < 0.05) allelic imbalances, with an allelic expression ratio of 2.82 observed in <it>BCL2L13 </it>(P = 0.004), 2.09 at <it>BID </it>(P = 0.001) and the most extreme case being at <it>MICAL3</it>, with an allelic expression ratio of 5.47 (P = 0.001). However, there was no correlation observed between the pattern of allelic expression and the genotype at rs2277831.</p> <p>Conclusions</p> <p>In the tissues that we have studied, our data do not support our hypothesis that the association between rs2277831 and OA is due to the effect this SNP has on <it>MICAL3, BCL2L13 </it>or <it>BID </it>gene expression. Instead, our data point towards other functional effects accounting for the OA associated signal.</p

Identification of new susceptibility loci for osteoarthritis (arcOGEN):a genome-wide association study

To access publisher's full text version of this article. Please click on the hyperlink in Additional Links field.Osteoarthritis is the most common form of arthritis worldwide and is a major cause of pain and disability in elderly people. The health economic burden of osteoarthritis is increasing commensurate with obesity prevalence and longevity. Osteoarthritis has a strong genetic component but the success of previous genetic studies has been restricted due to insufficient sample sizes and phenotype heterogeneity. We undertook a large genome-wide association study (GWAS) in 7410 unrelated and retrospectively and prospectively selected patients with severe osteoarthritis in the arcOGEN study, 80% of whom had undergone total joint replacement, and 11,009 unrelated controls from the UK. We replicated the most promising signals in an independent set of up to 7473 cases and 42,938 controls, from studies in Iceland, Estonia, the Netherlands, and the UK. All patients and controls were of European descent. We identified five genome-wide significant loci (binomial test p≤5·0×10(-8)) for association with osteoarthritis and three loci just below this threshold. The strongest association was on chromosome 3 with rs6976 (odds ratio 1·12 [95% CI 1·08-1·16]; p=7·24×10(-11)), which is in perfect linkage disequilibrium with rs11177. This SNP encodes a missense polymorphism within the nucleostemin-encoding gene GNL3. Levels of nucleostemin were raised in chondrocytes from patients with osteoarthritis in functional studies. Other significant loci were on chromosome 9 close to ASTN2, chromosome 6 between FILIP1 and SENP6, chromosome 12 close to KLHDC5 and PTHLH, and in another region of chromosome 12 close to CHST11. One of the signals close to genome-wide significance was within the FTO gene, which is involved in regulation of bodyweight-a strong risk factor for osteoarthritis. All risk variants were common in frequency and exerted small effects. Our findings provide insight into the genetics of arthritis and identify new pathways that might be amenable to future therapeutic intervention.Arthritis Research UK 1803

Rare variation at the TNFAIP3 locus and susceptibility to rheumatoid arthritis

Genome-wide association studies (GWAS) conducted using commercial single nucleotide polymorphisms (SNP) arrays have proven to be a powerful tool for the detection of common disease susceptibility variants. However, their utility for the detection of lower frequency variants is yet to be practically investigated. Here we describe the application of a rare variant collapsing method to a large genome-wide SNP dataset, the Wellcome Trust Case Control Consortium rheumatoid arthritis (RA) GWAS. We partitioned the data into gene-centric bins and collapsed genotypes of low frequency variants (defined here as MAF ≤0.05) into a single count coupled with univariate analysis. We then prioritised gene regions for further investigation in an independent cohort of 3,355 cases and 2,427 controls based on rare variant signal p value and prior evidence to support involvement in RA. A total of 14,536 gene bins were investigated in the primary analysis and signals mapping to the TNFAIP3 and chr17q24 loci were selected for further investigation. We detected replicating association to low frequency variants in the TNFAIP3 gene (combined p = 6.6 × 10−6). Even though rare variants are not well-represented and can be difficult to genotype in GWAS, our study supports the application of low frequency variant collapsing methods to genome-wide SNP datasets as a means of exploiting data that are routinely ignored

No evidence of an association between mitochondrial DNA variants and osteoarthritis in 7393 cases and 5122 controls.

OBJECTIVES: Osteoarthritis (OA) has a complex aetiology with a strong genetic component. Genome-wide association studies implicate several nuclear genes in the aetiology, but a major component of the heritability has yet to be defined at the molecular level. Initial studies implicate maternally inherited variants of mitochondrial DNA (mtDNA) in subgroups of patients with OA based on gender and specific joint involvement, but these findings have not been replicated. METHODS: The authors studied 138 maternally inherited mtDNA variants genotyped in a two cohort genetic association study across a total of 7393 OA cases from the arcOGEN consortium and 5122 controls genotyped in the Wellcome Trust Case Control consortium 2 study. RESULTS: Following data quality control we examined 48 mtDNA variants that were common in cohort 1 and cohort 2, and found no association with OA. None of the phenotypic subgroups previously associated with mtDNA haplogroups were associated in this study. CONCLUSIONS: We were not able to replicate previously published findings in the largest mtDNA association study to date. The evidence linking OA to mtDNA is not compelling at present

Evaluation of the genetic overlap between osteoarthritis with body mass index and height using genome-wide association scan data.

OBJECTIVES: Obesity as measured by body mass index (BMI) is one of the major risk factors for osteoarthritis. In addition, genetic overlap has been reported between osteoarthritis and normal adult height variation. We investigated whether this relationship is due to a shared genetic aetiology on a genome-wide scale. METHODS: We compared genetic association summary statistics (effect size, p value) for BMI and height from the GIANT consortium genome-wide association study (GWAS) with genetic association summary statistics from the arcOGEN consortium osteoarthritis GWAS. Significance was evaluated by permutation. Replication of osteoarthritis association of the highlighted signals was investigated in an independent dataset. Phenotypic information of height and BMI was accounted for in a separate analysis using osteoarthritis-free controls. RESULTS: We found significant overlap between osteoarthritis and height (p=3.3×10(-5) for signals with p≤0.05) when the GIANT and arcOGEN GWAS were compared. For signals with p≤0.001 we found 17 shared signals between osteoarthritis and height and four between osteoarthritis and BMI. However, only one of the height or BMI signals that had shown evidence of association with osteoarthritis in the arcOGEN GWAS was also associated with osteoarthritis in the independent dataset: rs12149832, within the FTO gene (combined p=2.3×10(-5)). As expected, this signal was attenuated when we adjusted for BMI. CONCLUSIONS: We found a significant excess of shared signals between both osteoarthritis and height and osteoarthritis and BMI, suggestive of a common genetic aetiology. However, only one signal showed association with osteoarthritis when followed up in a new dataset

Whole-Genome Sequencing Coupled to Imputation Discovers Genetic Signals for Anthropometric Traits

Deep sequence-based imputation can enhance the discovery power of genome-wide association studies by assessing previously unexplored variation across the common-and low-frequency spectra. We applied a hybrid whole-genome sequencing (WGS) and deep imputation approach to examine the broader allelic architecture of 12 anthropometric traits associated with height, body mass, and fat distribution in up to 267,616 individuals. We report 106 genome-wide significant signals that have not been previously identified, including 9 low-frequency variants pointing to functional candidates. Of the 106 signals, 6 are in genomic regions that have not been implicated with related traits before, 28 are independent signals at previously reported regions, and 72 represent previously reported signals for a different anthropometric trait. 71% of signals reside within genes and fine mapping resolves 23 signals to one or two likely causal variants. We confirm genetic overlap between human monogenic and polygenic anthropometric traits and find signal enrichment in cis expression QTLs in relevant tissues. Our results highlight the potential of WGS strategies to enhance biologically relevant discoveries across the frequency spectrum.Peer reviewe

The genetic epidemiology of joint shape and the development of osteoarthritis

Congruent, low-friction relative movement between the articulating elements of a synovial joint is an essential pre-requisite for sustained, efficient, function. Where disorders of joint formation or maintenance exist, mechanical overloading and osteoarthritis (OA) follow. The heritable component of OA accounts for ~ 50% of susceptible risk. Although almost 100 genetic risk loci for OA have now been identified, and the epidemiological relationship between joint development, joint shape and osteoarthritis is well established, we still have only a limited understanding of the contribution that genetic variation makes to joint shape and how this modulates OA risk. In this article, a brief overview of synovial joint development and its genetic regulation is followed by a review of current knowledge on the genetic epidemiology of established joint shape disorders and common shape variation. A summary of current genetic epidemiology of OA is also given, together with current evidence on the genetic overlap between shape variation and OA. Finally, the established genetic risk loci for both joint shape and osteoarthritis are discussed

Whole exome sequencing in an isolated population from the Dalmatian island of Vis

We have whole-exome sequenced 176 individuals from the isolated population of the island of Vis in Croatia in order to describe exonic variation architecture. We found 290 577 single nucleotide variants (SNVs), 65% of which are singletons, low frequency or rare variants. A total of 25 430 (9%) SNVs are novel, previously not catalogued in NHLBI GO Exome Sequencing Project, UK10K-Generation Scotland, 1000Genomes Project, ExAC or NCBI Reference Assembly dbSNP. The majority of these variants (76%) are singletons. Comparable to data obtained from UK10K-Generation Scotland that were sequenced and analysed using the same protocols, we detected an enrichment of potentially damaging variants (non-synonymous and loss-of-function) in the low frequency and common variant categories. On average 115 (range 93–140) genotypes with loss-of-function variants, 23 (15–34) of which were homozygous, were identified per person. The landscape of loss-of-function variants across an exome revealed that variants mainly accumulated in genes on the xenobiotic-related pathways, of which majority coded for enzymes. The frequency of loss-of-function variants was additionally increased in Vis runs of homozygosity regions where variants mainly affected signalling pathways. This work confirms the isolate status of Vis population by means of whole-exome sequence and reveals the pattern of loss-of-function mutations, which resembles the trails of adaptive evolution that were found in other species. By cataloguing the exomic variants and describing the allelic structure of the Vis population, this study will serve as a valuable resource for future genetic studies of human diseases, population genetics and evolution in this population